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Abstract 

It is an old speculation in physics that, once the gravitational field is success- 
fully quantized, it should serve as the natural regulator of infrared and ultraviolet 
singularities that plague quantum field theories in a background metric. 

We demonstrate that this idea is implemented in a precise sense within the 
framework of four-dimensional canonical Lorentzian quantum gravity in the contin- 
uum. 

Specifically, we show that the Hamiltonian of the standard model supports a 
representation in which finite linear combinations of Wilson loop functionals around 
closed loops, as well as along open lines with fermionic and Higgs field insertions at 
the end points are densely defined operators. 

This Hamiltonian, surprisingly, does not suffer from any singularities, it is com- 
pletely finite without renormalization. This property is shared by string theory. In 
contrast to string theory, however, we are dealing with a particular phase of the 
standard model coupled to gravity which is entirely non-perturbatively defined and 
second quantized. 



1 Introduction 

It is an old idea in field theory that once the gravitational field has been successfully 
quantized then it should serve as the natural regulator of matter quantum field theories. 
The argument is roughly that, since there is a fundamental length scale, namely the 
Planck scale £ p = y/hn where k is Newton's constant, the gravitational field should serve 
as an ultra-violet cut-off. The intuition coming from classical general relativity is that 
an elementary excitation of the fields whose energy exceeds the Planck mass will have an 
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energy density that is high enough in order for the excitation to become a black hole. Such 
tiny black holes, however, should evaporate within a Planck time scale into excitations of 
lower energy. The conclusion is that elementary excitations will have a 4-volume which 
is larger than or equal to 

In the present article we show a precise realization of this idea within the framework 
of canonical quantization of matter-coupled four-dimensional Lorentzian quantum gravity 
in the continuum. In brief, it is actually possible to find a representation in which all 
matter Hamiltonians become, after suitable regularization, densely defined operators. 
Upon removing the regulator no renormalization is necessary and so the theory is, just 
like string theory, completely finite. In contrast to string theory, however, our approach 
is fully non-perturbative and starts from the second quantized field theory. In particular, 
our framework is mathematically completely rigorous, we have a well-defined Hilbert 
space and all the matter Hamiltonian operators are densely defined on it. Approximation 
schemes, if necessary to solve the theory, would not be formal perturbation series with little 
if not no control on the error, but approximation schemes with full control of convergence 
issues just like in usual quantum mechanics. 

The intuitive picture that arises from the Hilbert space we choose is as follows : 
The elementary excitations of the gravitational and gauge fields are concentrated along 
open or closed strings while those of the fermion and scalar fields are located in the 
endpoints of the open strings. The Hamiltonian of the standard model and gravity, in 
this diffeomorphism invariant phase of the full theory, act by creating and annihilating 
those excitations which reminds of a non-linear Fock representation. 

It should be stressed from the outset, however, that the string enters here as a com- 
pletely kinematical object and unlike in string theory does not acquire any dynamical 
properties. It is just a label for the state, in fact, the same label that one uses for the 
familiar Wilson loop functionals that one knows from lattice gauge theory. Moreover, in 
contrast to string theory, the strings that acquire physical importance have necessarily 
singularities, i.e. they intersect in an arbitrarily complicated, non-differentiable, manner. 



The plan of the paper is as follows : 

In section 2 we recall the quantum kinematics of the canonical approach from |Q 0, 
|], |^, [J for the gravitational and gauge sector and from [|7| for the Fermion and Higgs 
sector (see also || |9| for earlier work on the Fermion sector which, however, is described 
by a Hilbert space with an inappropriate inner product). 

In section 3 we come to the regularization of the matter Hamiltonians which is very 
similar in nature to the one performed for the gravitational field, in fact the techniques 



used extend those introduced in fll|, |12|. Roughly, what we do is to introduce an ultraviolet 
cut-off by triangulating the spacelike hypersurfaces £ of the four-dimensional spacetime 
M = R x E and then to take the continuum limit. We show that it exists and are able 
to precisely display the action of the continnum operator. At no stage we encounter any 
singularities, these final operators do not require any renormalization. 

In the first subsection we regulate the QCD Hamiltonian for any compact gauge group 
G. Not surprisingly, the electric part of the Yang-Mills Hamiltonian becomes a sum of 
Laplace-Beltrami operators on G and therefore is sensitive to the colour of the state while 
the magnetic part creates and annihilates new excitations, that is, it creates new Wilson 
loop functions. 

In the next subsection we address the fermionic term. The fermionic Hamiltonian oper- 
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ator removes fermionic excitations from open string endpoints and creates new excitations 
on the open string. Quite surprisingly, these fermionic excitations are very different from 
those discussed in || || the underlying reason being the faithful implementation of the 
reality conditions. 

In the following subsection we discuss the regularization of the Higgs Field Hamilto- 
nian. The action of that Hamiltonian is analogous to the one of the gauge field. 

Finally, in the last subsection we outline a general procedure for regulating a rather 
general class of Hamiltonians whose corresponding Hamiltonian density has a density 
weight of one. 

At this point the reader will be puzzled what role the gravitational field still plays. 
As will become obvious from the details, it enters the stage simultanously in two different 
ways (remember that the gravitational field couples to matter always through the three- 
metric q a b or the co-triad e l a of £ (in addition, fermions couple to the extrinsic curvature 
as well)) : 

• UV Singularities 

Recall that one may control the usual ultraviolet singularities in terms of point- 
splitting regularizations of operator-valued distributions multiplied at the same 
point. For instance, we may have a singular square of operator- valued distribu- 
tions of the form F(x) 2 which one may regulate by point splitting / d 3 yf e (x — 
y)/e 3 F(x)F(y) and \im t ^ f t (x) / e 3 = 5(x). Notice that we automatically have bro- 
ken diffeomorphism covariance because the points x, y are split by a background 
metric. It turns out that the point splitting volume e 3 is absorbed by a certain 
gravitational operator, built from q^ that measures the volume of spatial regions. 
This is intuitively reasonable because the three-dimensional coordinate volume e 3 
cannot be measured by a fixed background metric in a diffeomorphism invariant 
theory like general relativity but must be measured by the dynamical metric q a b 
itself ! This volume operator therefore must enter the final expression of all mat- 
ter Hamiltonians. The formalism itself predicts how it enters, we do not have to 
postulate this, of course, up to ordering ambiguities. Since this volume operator 
turns out to be densely defined on the Hilbert space the UV singularity 1/e 3 is re- 
moved by coupling quantum gravity, without renormalization, thereby reinstalling 
diffeomorphism covariance. 

• IR Singularities : 

This volume operator turns out to have a quite local action, it vanishes everywhere 
except at points where the string ends or starts. This is also an unexpected pre- 
diction of the formalism. It is this feature which makes the Hamiltonian operators 
densely defined without that we have to introduce an infra-red (infinite volume) 
cut-off. 

In a sense, it is the volume operator which is the natural regulator of the matter quantum 
field theories by serving as a dynamical ultra-violet and infra-red cut-off ! 

And in accordance with what we said at the beginning of this section, the volume is 
quantized with discrete spectrum, the quantum of volume being indeed of order £ 3 (see 
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21 1) ! In particular, it will become obvious in the course of the construction that 
matter field theories are ultraviolet and infrared 

a) divergent without gravity and need to be renormalized but are 

b) convergent with gravity without that renormalization is necessary. 

In section 4 we perform various consistency checks on the theory, for instance, that we 
do not encounter quantum anomalies when computing commutators. This can be done 
only by restricting to the diffeomorphism invariant subspace of the Hilbert space in which 
we are interested only. Expectedly, it is also in this context only that we can remove yet 
another ultra-violet regulator for the connection field which enters in terms of a triangula- 
tion of E. When refining the triangulation ad infinitum we find that the continuum limit 
exists and yields well-defined operators on this diffeomorphism invariant Hilbert space. 

We also address the question whether the operators obtained are positive semi-definite, 
at least on the kinematical Hilbert space which, in view of some kind of "quantum dom- 
inant energy condition", would be a re-assuring result because the matter Hamiltonian 
constraint plays the role of the timelike-timelike component of the energy momentum 
tensor. We find that, for each matter species separately, the answer is regularization 
and factor-ordering dependent. We clarify the meaning of this result and point out that 
what only is important is that the total (ADM) energy is non-negative (see [jllj for the 



pure gravity case). Finally, we comment on the general construction of solutions to the 
full Hamiltonian constraint and demonstrate non-triviality of the theory by displaying an 
uncountably infinite number of rigorous simple solutions. 

In appendix A we derive the Dirac-Einstein canonical action in manifestly real form 
and in terms of the real connection variables that have proved successful in quantizing 
the source-free gravitational field in [EJ . This, to the best of our knowledge, has not yet 
been done in the literature. 



2 Preliminaries 

We begin by describing the field content of the theory. 

The topology of the four-dimensional manifold is chosen, as always in the canonical 
approach, to be M = R x E where E is a smooth 3-manifold which admits smooth 
Riemannian metrics. 

On E, there is defined a co-triad field e\ where a,b,c,.. denote tensorial indices and 
i,j, k, .. denote su(2) indices. From this field the 3-metric is derived via q a b = e l a e\. Fur- 
thermore, we have a field K\ from which the extrinsic curvature K a f, of E is derived via 
K a b = sgn(det((ej))).?Qe&. It turns out that the pair (K l a , E?/k) is a canonical one on the 
gravitational phase space where Ef := ^e abc eijke J b e^. 

The Hamiltonian constraint (or Wheeler-DeWitt constraint) of general relativity takes a 
quite complicated form in terms of these variables, thus it was natural to assume that 
for purposes of quantizing canonical gravity it is mandatory to cast the theory into poly- 
nomial form. The famous discovery due to Ashtekar [[03] is that this indeed possible by 
performing a certain canonical transformation on the gravitational phase space. However, 
this transformation comes at two prizes : 

1) The Hamiltonian constraint is polynomial only after rescaling it by Jdet((q a b))- This 
is bad because the constraint now adopts a density weight of two which rules out a dif- 
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feomorphism covariant regularization and will require a troubelsome multiplicative renor- 
malization. This will become apparent in section 3. 

2) The so-called canonical Ashtekar variables i^A l a = T l a — iK l a , iEf/n), where T l a is the 
spin-connection of e*, are complex-valued. This is bad because the Ashtekar connection 
^A^ is the connection of a principal SL(2, C) bundle, that is, the gauge group is non- 
compact and makes the rich arsenal of techniques that have been developed for gauge 
theories with compact gauge group inaccessible. 

There have been two quite different proposals to deal with problem 2). First of all, in [|14[ 
it was pointed out that one can also use a real-valued Ashtekar connection A\ = T l a + K l a 
at the prize of living with a fairly complicated Lorentzian Hamiltonian constraint. The 
virtue is that this at least opens access to the techniques developed in [[[], ||, [| £|, [|, || and 
equips us with a Hilbert space structure that faithfully implements the reality conditions. 
Restricted to Euclidean gravity this was also proposed in fll5] . 



The second proposal is to perform a Wick rotation on the canonical phase space ||16|| . The 
virtue of this is that one can start by quantizing the Euclidean Hamiltonian constraint in 
terms of the real canonical variables {A^E^/k) which takes care of the reality structure 
of the theory and keeps the constraint polynomial. In a second step then one would per- 
form a Bargman-Segal kind of transform to the Lorentzian theory described by complex 



valued connections (compare also a modified procedure jL?| which could enable one to 



stay purely within a real connection theory). The drawback is that the generator of the 
Wick transform adopts a quite complicated form which made it hard to imagine how one 
would be able to quantize it (see, however, |12| for a proposal for a self-adjoint operator). 



Apart from the problems mentioned, both proposals still suffer from the problem 1) de- 
scribed above. 

In ]TT|, |T2j| a novel technique was introduced which solves both problems 1),2) in one 



stroke and on top defines the generator of the Wick rotation transform. The resulting 
Lorentzian Hamiltonian constraint is densely defined, anomaly-free and one has a Hilbert 
space that incorporates the correct reality conditions. 

This paper is devoted to the extension of this technique to the non-gravitational sector. 

Let G be an arbitrary compact gauge group, for instance the gauge group of the stan- 
dard model. Denote by J, J, K, ... Lie(G) indices. We introduce classical Grassman- valued 
spinor fields r\ = (t/a^) where A,B,C, .. denote indices associated with the gravitational 
SU(2) and /z, z/, p, .. with the group G. The fermion species 7] transforms like a scalar and 
according to an irreducible representation of SU (2) x G. It turns out that in its mani- 
festly real form (the associated conjugation is just complex conjugation for non-spinorial 
variables and for spinorial fields it involves a cyclic reversal of order in products) the 



most convenient description of the constraints is in terms of half-densities £ := y det(q)r). 

The momentum conjugate to £a,h is then just given by tca,^ = i^A,u an d the real-valued 
gravitational connection is given by A l a = T l a + K l a just as in the source free case. As 
we will see in appendix A, the connection is real only if we use the quantities £ with 
density weight 1/2, if we would use the scalar variables 77 as in [[| then the gravitational 



connection becomes by the argument given in [18 



4 = r* B + Kl + — = e fc (2.1) 
4Jdet(g) 



which is complex valued and therefore makes the techniques in || ||, f|, ^ || inaccessable. 
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Notice that it is no lack of generality to restrict ourselves to just one kind of helicity : If 
we have several fermion species then we can always perform the canonical transformation 
~^ (^£5 e where e is the spinor-metric, the totally skew symbol in two dimensions. 
Notice that there is no minus sign missing because we take the fermion fields to be anti- 



commuting, the action is form- invariant under this transformation [18 . 

In the gauge sector we have canonical pairs (A^,E_j/Q 2 ) where the first entry is a 
G connection and the second entry is the associated electric field, Q is the Yang-Mills 
coupling constant. Finally, we may have scalar Higgs fields described by a canonical 
pair ((pijp 1 ) transforming according to the adjoint representation of G. Without loss of 
generality we can take these as real valued by suitably raising the number of Higgs families. 
Here and in what follows we assume that indices /, J,K,.. are raised and lowered with 
the Cartan-Killing metric 5u of G which we take to be semi-simple up to factors of U(l). 

We could also introduce Rarity-Schwinger fields and make everything supersymmetric 
but since this will not add new features as compared to the ordinary spinorial action, we 
refrain from doing so. 

With this field content we then have the following Lorentzian Hamiltonian constraints : 

1 



H Einstein = tv(2{[K a , K b ] - F ab }[E\ E b }) + X^det(q) 

KWdet(g) 



H Di rac = Ef 1 \m T T{D a i + V^ni) + ^A> T £ + c.c] 



2Qydet( 9 ) 



H Higgs = -{ y + y/det(q)[q ab (V a <f> I )(V b <f> I )/K + P^j)/^ 2 )]). (2.2) 
z nJdet(q) 



Here we have denoted by 7j the generators of the Lie algebra of su(2) with the conven- 
tion [ri,Tj] = €ijkT k , F ab is the curvature of A a (one can check that all the constraints 
remain form-invariant under the canonical transformation that turns the fermions into 
half-densities, see appendix A), T> is the covariant derivative with respect to SU{2) x G, 
that is, with respect to w a := A a + A a and B_ a is the magnetic field of the Yang-Mills 
connection. We have included a cosmological constant (A) and P denotes an arbitrarily 
chosen gauge invariant function of the Higgs field (not including spatial derivatives), the 
Higgs potential. Notice that we have rescaled the Higgs field by yH in order to make it 
dimensionless. 

The unfamiliar terms in the Dirac Hamiltonian proportional to the total derivative and 
K l a arise because 1) we are dealing with half densities rather than scalars and 2) we couple 
the real connection A to the spinor fields while in the traditional approach it is naturally 
the complex valued (anti-)self-dual part of the spatial projection of the spin-connection 
that couples to them. Thus, these additional terms are the required correction terms if 
we describe the theory in the variables we chose. The interested reader is referred to 
appendix A in order to see how these corrections come about. As usual, the "c.c." means 
involution (complex conjugation for complex valued fields and an additional reversal of 
order is implied for the Grassman valued fields). 

In ( |2.2| ) we have written only one family member of the possibly arbitrary large family of 
field species, in particular, we can have an arbitrary number of gauge fields all associated 
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with different gauge groups and associated "quarks" and "Higgs" fields and transforming 
under different irreducible representations of 577(2) x G. However, we will not deal with 
these straightforward generalizations and consider only one species of fermions or Higgs 
fields respectively which transform under the fundamental representation of both SU(2) 
and G or the adjoint representation of G respectively. Also, one could easily deal with 
a more complicated "unified gauge group" which is not of the product type SU (2) x G 
but contains it as a subgroup. However, for simplicity and because one does not expect 
a unification of the gauge group of the standard model and the gauge group underlying 
the frame bundle, we refrain also from treating this more general case. 
This furnishes the description of the classical field content. 

We now come to the quantum theory. We can immediately apply the techniques of 
|l|, 0, 0, f|, H |6| to write down a kinematical inner product for the gravitational and 
Yang-Mills sector that faithfully incorporates all the reality conditions. We get a Hilbert 
space L 2 (Asu(2) x Ac,dfiAL,su(2) <8> g^al.g) where the index "AL" stands for Ashtekar- 
Lewandowski measure and the group index indicates to which gauge group the Ashtekar- 
Lewandowski measure is assigned. The reader interested in the constructions and tech- 
niques around the space of generalized connections modulo gauge transformations is urged 
to consult the papers listed. In particular, the probability measure Hal is very natural 
and diffeomorphism invariant. If we are interested in gauge invariant functions of connec- 
tions alone, then the space of generalized connections A can be replaced by the space of 
generalized connections modulo gauge transformations A/Q. 

The extension of the framework to Higgs and fermionic fields is not entirely straight- 
forward : 

Let us first focus on the Higgs field. Assume that we choose (f>i(x) as our basic configura- 
tion field variable. As argued in , in a diffeomorphism invariant theory this assumption 
has consequences which leads to inconsistencies. Basically, the problem is the following : 
The variables <pi{x) are real- valued and thus there does not exist a translation invariant 
measure on the space of these 0/'s. For a quantum field theory in a fixed background 
there is no problem, a natural kinematical measure that incorporates the reality condi- 
tions is a Gaussian measure leading to a usual Fock Hilbert space. However, a Gaussian 
measure for a scalar field, rigorously defined through its covariance, is always background 
dependent or, in other words, cannot be diffeomorphism invariant ! An intuitive way to 
see this is by recalling that the covariance is determined by the characteristic functional 
exp(J d 3 x J d 3 yCu(x, y)f I (x)f J (y)) of the measure which in turn is the expectation value 
of exp(i / d 3 xf T (x)(f)i(x)) where f 1 are some test functions. However, the fact that 4>i is 
a scalar implies that the kernel Cjj(x,y) of the covariance is a density of weight one and 
therefore the characteristic functional is background dependent. See |7j for more details. 

Thus we need a new approach which does not use a Gaussian measure and therefore 
we must not use <pi as a basic variable but some variable that is valued in a bounded 
set. This motivates to use the variables U(v) := exp(0/(t>)r/) quite in analogy with the 
holonomy for a connection and we will call them "point holonomies" . Point holonomies 
are G- valued and, since G is compact, its matrix elements are therefore bounded. In |7| 
we construct a representation in which the U(v) are promoted to unitary operators (since 
we can replace G by a unitary group by the theorem due to Weyl that any compact group 
is equivalent to a unitary one). If we are dealing not with a Higgs field but just with a real 
scalar field then we may use U(v) = e^W. The Hilbert space to be used is surprisingly 
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simple to describe : there is a certain space U of generalized Higgs fields which turns 
out to be in bijection with Fun(£, G), the space of all functions from £ to the gauge 
group. That is, a typical such function is a "wild", arbitrarily discontinuous function, it 
is a wild Higgs field. On that space we have a measure /ij/ which is a rigorously defined 
cr-additive probability measure on U which is formally given by the uncountable direct 
product dfijj^cf)) := Y[ v &t, d^H(U(v)) where fin denotes the Haar measure on G. The 
Hilbert space is then the corresponding L 2 (W, d^u) space and one can show that this 
is the unique Hilbert space selected by the adjointness relations, once we have chosen the 
space U as the quantum configuration space. Expectedly, the mathematical description 
is very similar to the one for gauge fields |T], |2], |3|, £|, |5|, |6| . 

Next we come to the fermion fields £ which, as explained above, have density weight 
1/2. 

It turns out |7j] that the faithful implementation of the reality conditions forces us to work 
in a representation in which the objects 

become densely defined multiplication operators. Here x<l{ x i u) is t ne characteristic func- 
tion of a box of Lebesgue measure e 3 and center x. The 9 are by inspection scalar 
Grassman-valued quantities because the 5 distribution is a density of weight one. In cal- 
culations it is understood that the e — > limit is performed only after the manipulation 
under consideration is performed [[/J. 

Consider then the n = 2d Grassman variables 9i(v), A = 1, 2, fj, — 1, .., d where d denotes 
the dimension of the fundamental representation of G. Here we have have introduced a 
compound symbol i instead of Afi to simplify the notation. These variables coordina- 
tize together with their conjugates the superspace S v at point v. Since Grassman fields 
anti-commute, any product of more than n of these Oi(v), i = l,..,n will vanish. The 
vector space of monomials of order k is n\/(k\(n — A;)!) dimensional where k — 0,1, ..,n 
and the full vector space Q v built from all monomials has dimension 2 n . The quantum 
configuration space is the uncountable direct product ("superspace") S := Il^gs &v an d 
in order to define an inner product on S it turns out to be sufficient to define an inner 
product on each S v coming from a probability "measure". The "measure" on 5*,, is a 
modified form of the Berezin symbolic integral [19| : 

dm(9,9) = dBdQe ee and dm v = ^ =1 dm(6 l (v),e i (v)). 

The fermionic Hilbert space is then simply given by 

H F = L 2 (S, dfi F ) = ® V (z S L 2 (S v , dm v ) 

where "F" stands for fermionic and it is understood that we integrate only linear combi- 
narions of functions of the form fg where /, g are both holomorphic (that is, a function on 
S which depends on 6i{v) only but not on 6i{v)). As a result, the integral of any function 
of the type /*/, where / is any holomorphic function, is strictly positive and so we have 
an inner product. This inner product, when restricted to one point v, is easily seen to be 
the standard inner product on Q v when viewed as the vector space of exteriour forms of 
maximal degree n. Thus, Tip is a space of holomorphic square integrable functions on S 
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with respect to djj, F . The Fermion measure /ip is easily seen to be gauge and diffeomor- 
phism invariant. 

The reader is referred to |7[ for a more complete treatment where it is also shown that the 
reality condition = — m is faithfully implemented in the inner product. The reader will 
find there also an extension of the framework to the diffeomorphism invariant subspace 
of the Hilbert space. 

Let us summarize : the Hilbert space of (not necessarily gauge invariant) functions of 
gravitational, gauge, spinor and Higgs fields is given by 

H := L 2 (A SU(2) , d^ AL (SU (2))) ® L 2 (A G , dfi AL (G)) ® L 2 (S, d^ F ) ® L 2 (U, d^u). 

The Hilbert space of gauge invariant functions will be just the restriction of Ti to gauge 
invariant functions. It turns out that, because our total measure is a probability measure, 
gauge invariant functions will be still integrable with respect to it, in other words, "the 
gauge group volume" equals unity in our case ! 

A natural gauge invariant object associated with spinor fields, Higgs fields and gauge 
fields are "spin-colour- network states" J7|. By this we mean the following : Let 7 be a 
piecewise analytic graph with edges e and vertices v which is not necessarily connected 
or closed. By suitably subdividing edges into two halves we can assume that all edges are 
outgoing at a vertex. Given a (generalized) connection u a = A a + we can compute the 
holonomies h e (A), /i e (A), H e {uj) = /i e (^4)/leGl)- With each edge e we associate a spin j e 
and a colour c e corresponding to irreducible representations of SU (2) and of G respectively 
(for instance for G = SU (N), c e is an array of iV — 1 not increasing integers corresponding 
to the frame of a Young diagramme). Furthermore, with each vertex v G Vi^y) we associate 
an integer n v , yet another colour C v and two projectors p v , q v . Here ^(7) denotes the 
set of vertices of 7. The integer n v corresponds to the subvector space of Q v spanned by 
monomials of degree n v . 

Likewise, the colour C v stands for an irreducible representation of G, evaluated at the 
point holonomy U(v). The projector p v is a certain SU(2) invariant matrix which projects 
onto one of the linearly independent trivial representations contained in the decomposition 
into irreducibles of the tensor product consisting of 

a) the n v — fold tensor product of fundamental representations of SU(2) associated with 
the subvector space of Q v spanned by the monomials of degree n v and 

b) the tensor product of the irreducible representations j e of 577(2) of spin j e where e 
runs through the subset of edges of 7 which start at v. 

Likewise, the projector q v , repeats the same procedure just that S77(2) is being replaced 
by G and that we need to consider in addition the adjoint representation associated with 
C v coming from the Higgs field at v. Now we simply contract all the indices of the tensor 
product of 

1) the irreducible representations evaluated at the holonomy of the given connection, 

2) the fundamental representations evaluated at the given spinor field and 

3) the adjoint representations evaluated at the given scalar field, 

all associated with the same vertex v, with the projectors p v ,q v in the obvious way and 
for all v G ^(7). The result is a gauge invariant state 
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which we will call a spin-colour-network states because they extend the definition of the 
pure spin-network states which arise in the source- free case (e.g. ||). 
These spin-colour-networks turn out to be a basis for the subspace of gauge invariant func- 
tions. They are not orthonormal, but almost : we just need to decompose the fermionic 
dependence into an orthonormal basis for each of the Q v Q|. 

This furnishes the summary of the quantum kinematics. We now turn to the quantum 
dynamics. 



3 Regular izat ion 

The regularization of the Wheeler-DeWitt Hamiltonian constraint was carried out in 



1% |12 |. We therefore can focus on the remaining Hamiltonians. 



3.1 Gauge sector 

We begin by looking at the electric piece Q 2 H§ M (N) = J d 3 xN — ^ ^ E^E b j. Recall from 



11, 12] that the following identity was key 



V ^ = JEf = 2sgn(det((4))K (3.1) 

where V = J d 3 xyjdet(q) is the total volume of the hypersurface (in the asymptotically 
flat case the appropriate definition of the functional derivative of V involves a certain 



limiting procedure |TT 



Let now e be a small number and let Xe( x , v) — Ila=i 0(e/2 — \x a — y a \) be the characteristic 
function of a cube of coordinate volume e 3 with center x. That is, we have chosen some 
frame andd therefore broken diffeomorphism covariance in the regularization step. We are 
going to remove the regulator later again and also recover diffeomorphism covariance. Also 
let V(x, e) := J d 3 yXe(%, 2/)\/det(g) be the volume of that box as measured by q a b- Then, 



since lim e ^ jsXe(x,y) = $(x,y) we have lim^o jsV(x, e) = ^det(q)(x). It is also easy to 
see that for each e > we have that 6V/5E°-(x) = SV(x, e)/5E^(x). The simple trick that 
we are going to use is as follows : let /, g be some integrable functions on £ with respect to 
Lebesgue measure. Thus f(x) = lim e ^ J d 3 yf(y)x e (x,y) =: lim e ^ \f{x, e) and sim- 
ilar for g. Then lim e ^ [/(a;, e)/sf(|/, e)] = \\m e ^ [{ f(x, e)/e 3 }/{g(y, e)/e 3 }] = f(x)/g(y), 
that is, the two singular factors of e 3 cancel each other in the quotient. 
With this preparation it follows that we have the following classical identity 

2k 2 Q 2 H* m (N) 

= limi f d s xN{x) {^x)V} I d 3 ( } { ^ y ) V } E](y) 



2^det(g)(x) J 2^det(g)(2/) 

Ami /^mw^ 

e ^ 0e ^ 2^det(g)(x) J 2tfdeb(q)(y) 
hm- / d xN(x) Ej(x) / d yXe{x,y) . E^y) 
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lim / d 3 xjv(x){4(x), v/k(x,6)}^(x) / Ax e (^,y){4(y),\A%^)}^/(y) 



(3.2) 



which demonstrates that we can neatly absorb the annoying 1 / ^/detq into a Poisson 
bracket, of course at the prize of breaking gauge invariance at finite e. The removal of 
the divergent factor 1/e 3 has occured precisely because we kept the density weight of the 
constraint to be one ! 
We have similarily for the magnetic term 

2Q 2 k 2 H* m (N) = lim j rf 3 xiY(x){^(x) )V /V(x, e )}^(x) J d 3 y X e(x,y) x 



x {Al(y),^V(y,e)}B>(y). (3.3) 

We come now to the quantization of the Yang-Mills Hamiltonian constraint. This will 
be somewhat different for the electric and the magnetic part so that we describe them 
separately. Notice that we have no factor ordering problem at all as far as the question, 
whether to order the gravitational or the gauge theory variables to the left or to the right, 
is concerned. 

Let us then start with the electric part. We choose to order the Yang Mills electric fields 
to the right and replace E_} — > —ih5/5A^, V — > V (V(R), for an arbitrary region R, 
was described in [p0| , pT| 22, |23|]) and Poisson brackets by commutators times l/ih. If 



yhK^rrip = tJH/k denote Planck length and mass respectively then we obtain on a 
function / cylindrical with respect to a graph 7 the following result (ccq = %Q 2 is the 
dimensionless fine structure constant )Q 

-H?&(N)f 



Y^^L J ^xN(x)[^ a (x),y/v(x,e)] J d 3 y X e(x,y)[Ai(y)Jv(y,e)] x 



e,e' ~"p 

„//j./\\ -lb/j.1 



x / dt5{x,e{t))e a {t) dt'5{y,e'{t'))e' b {t') x 
Jo Jo 

x {[tr(/^(0,t)T//t e (t, l)d/dh e (0, l))tr(/i e ,(0, £>/&e'0'> 1)0/0M°> 1))] 
+5 ei e^(t' - t)tr(/i e (0,t)r / /i e (t,t / )^ e (^, 1)9/^(0, 1)) 
+0(t - t / )tr(4e(0,t / )rA(t / ,t)r / /i e (t, 1)9/9^(0, 1))]}/ 

= E^Srf'dt t dt'N{e{t)) Xe {e{t),e'{t')) x 
ee , 2^ Jo 7o 

x K(e(t))e a (t), ^((e^.c)]^^))^^), V^W] x 
x {[tr(^(0, t)rj^(t, l)d/dtu(p, l))tr(h e ,(0,t')TjhAt', l)d/dh e ,(0, 1))] 
+5 e , e >[0(t' - OtrCkCO^Tj/^t,*^^^, 1)9/9^,(0, 1)) 

- t , )tr(4 e (0,t / )^ e (t',t)rA(t, l)d/dhe(0, 1))]}/. (3.4) 



^ere and in the regularizations that follow we are going to apply the operator first only to functions 
of classical (i.e. smooth) fields and then extend the end result to the quantum configuration space 
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Here we have used the step functions 9(t) — 1 if t > and otherwise. The negative sign 
in (|3.4j) stems from the (— i) 2 coming from the two factors of the electrical field. 



The next step consists in replacing the integrals by Riemann sums, suggested by the 
expansion [h e (t,t + St),0] = 5te a (t)[A a (e(t)), O] + o(5t 2 ) for an arbitrary operator O. 
We choose an arbitrary partition of the interval [0, 1] into n intervals with endpoints 
t k , k = 0, ..,n which we can since the Riemann integral is independent of the partition 
that defines it (here we used that for the moment being we deal with smooth connections). 
Since the formula {A l a (y), V(x, e)} oc e l a (y) is always true provided that Xe( x ,y) = 1 we 
have that {h e (t k , tk+i), V(e(t k ), e)} ^ whatever partition we choose. We may therefore 
choose for given 7 the partition, given e, such that the following two conditions hold : 

1) min eifcja (|e a (t fc ) - e a (t fc _i)|) > e and 

2) min e ^ e / ifc+i> o, a (|e a (t fc ) - e ,a (t,)|) > e. 

Thus, the partition is as fine as we wish but only so fine that that (t k — t k+ i) /e is at least 
of order o(l). 

We then can replace (|3.4|) by 



-H?&{N)f 

e,e' Zl p k,l=l 

x ti(h e (tk, t k+1 )[h e (t k , tfc+i) -1 , \/v(e(t k ), e)] x 
x h e ,(t h t l+1 )lh e ,(t h t l+1 )- 1 , y/v(e(ti),e)]) x 

x {[tr(^(0, t k )rjh e (t k , l)d/dh e (0, l))tr(^ e ,(0,^)T/k'(*/, 1)9/3^(0, 1))] 
+5 e A0{k - t fe )tr(^ e (0,t fc )r / /i e (t fc) ^)T/A e (t i , 1)9/^(0, 1)) 
+6(t k - *,)tr(&e(0, tdnhefa, t k )Tjh e (t k , l)d/dh e (0, l))}}f. (3.5) 

Consider first the terms with e 7^ e'. Then for sufficiently small e we get Xe{ e (t k ), e'iU)) = 
unless e, e' intersect each other. We have set up the problem in such a way that they then 
must intersect in a vertex e(0) = e'(0) = v of the graph 7. Now, by condition 2) on our 
partition we also obtain that Xe( e (^fc), e '{ti)) — unless t k = ti = 0. 

If e = e' then Xe(e(t k ), e {U)) = unless k = I by condition 1) on the partition. Now we 
make use of the fact that we can commute the gravitational operators with the Yang-Mills 
operators. We then find out that 

tr(/i e (t fc ,t fc+ i)[/i e (t fc ,t fc+ i) _1 , y/v((e(t k ),e)] x 
x h e (t k , tk+i)[h e (tk, tk+iy 1 , \/v((e(t k ), e)])f 

= -tr([h e (t k ,t k+1 ), \/v((e(t k ),e)][h e (t k ,t k+1 )-\ \/v((e(t k ),e)])f (3.6) 

vanishes unless k = because the volume operator V(x, e) annihilates a state unless there 
is a vertex in the region corresponding to the e— box around x and because, by definition, 
only the starting point of an edge is a vertex of the graph. But by definition 6(0) = 0. 
We thus conclude that for sufficiently small e we obtain 



H' YM {N)f = £ £ tT(h e (0,e)[h e (0,e)-\^V(v,e)]x 

P v£V(-y) v£er\e' 



h e/ (0,e)[h e ,(0,e)-\^V(v,e)])2[lXl'f (3-7) 
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which does not depend on the details of the partition any longer because of which we could 
replace t% by e ! Here we have defined the right-invariant vector fields X J '(g) = tr(rjgd/ dg) 
and 2Le = 2L(he)- Notice that the final form of the operator fl3.7| ) is manifestly gauge 
invariant. Also, we could actually replace V(v, e) by V because the commutator [h e , V] 
equals [h e , V(e(0), R)} for any arbitrarily chosen neighbourhood of v — e(0), see |TT|, p^| . 
In particular, it equals [h e , V e (o)] where V v denotes the volume operator at a point v. This 
operator is defined on any cylindrical function / 7 as follows : consider an arbitrary finite 
contractable neighbourhood R of v, denote by R t any homotopy with R x = R, R = {v} 
and evaluate V(R t )f-y. By the properties of the volume operator, the vector V(R t )f 1 
is constant for all t < i 7 for some value t 7 > which depends on 7 whenever R t is so 
small that v is possibly the only vertex of 7 contained in R t and the vector is moreover 
independent of R and the homotopy. This vector is denoted by V vn f 1 . The family of 
operators {^,7)7 so defined is consistently defined because V(R) is and therefore qualifies 



as the cylindrical projection of an operator V v (see [27] for a different definition in terms 
of germs of analytical edges). 

This replacement of V(v, e) by V v is possible only in the quantum version where the square 
root of volume operator is defined via its spectral resolution and automatically takes care 
of its local action while in the classial computation we would have to keep the e. Thus the 



only e dependence of (|3.7| ) rests in the holonomies h e (0,e). But since the operator (|3.7|) 



is gauge invariant, by an argument given in [24| even the remaining e— dependence drops 



out as follows : we define for each edge e of the graph incident at v a segment s(e) also 
starting at v but not including the other endpoint of e. After evaluating the operator on 
a state, the dependence on s(e) automatically drops out. 

The final expression for the electrical part of the Yang- Mills Hamiltonian is then given by 
H? M (N)f = -^E E Nv E x 



X Hhs(e)[K^ VKl/MeO^e')' VK])2&£/ (3-8) 

where N v = N(v). 

Notice that we have exchanged the limits of taking e — >• and the limit of refining the 
partition ad infinitum. However, one could have arrived at (|3.8j) also differently : let 



5 := inffc(tfc — tfc_i). Make e in ( |3.2|) y-dependent, that is, e(y) = p(y)5 where p(y) > 
d 3 y almost everywhere and such that the conditions 1),2) on the partition hold (with e 
replaced by e(y)) at y = e(ife)). Then instead of taking e sufficiently small and the partion 
small but still finite we make e dependent on 5 in this sense and just take S sufficiently 
small but still keep it finite. The result ( |3.8| ) is the same by construction, just that we did 
not need to take any limits and so the questionable interchange of limiting procedures is 
unnecessary. Now, since ( |3.8j ) actually is independent of 5, no limit needs to be taken. 
On the other hand, this latter regularization scheme is, in contrast to the former scheme, 
state-dependent although the final operator is state-independent as we will see in the next 
section. 

Let us now turn to the magnetic part. In this case we need to introduce a triangula- 
tion of S just as in [11, |l2j in order to define its regularization. Taking over the no- 



tation from [[11], for the triangulation of E, for each vertex of 7 and each triple of 
edges e, e', e" we introduce tetrahedra A with basepoint v (A) = v and incident segments 
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Si (A), i = 1, 2, 3 where there is a one to one map between the segments s(e), s(e'), s(e") 
as defined above and the three Sj(A). We have assumed that 6 a (, c s" > 0. We will 
denote the arcs of A that connect the endpoints of Sj(A),Sj(A) by ay (A). Finally we 
have loops := s, o a^- o sj 1 . 

Given first of all any triangulation of E it is immediate to see that the magnetic part of 
the Yang-Mills Hamiltonian can be written, by the same trick that we used for the electric 
part, as follows 

H* M (N) = hrn-Li: f N(x){A\x), y/v(x, e)} A F_\x) x 

x / Xe (x,y){A i ( 2/ ), v /v^)}AF / ( 2/ ). (3.9) 



Now notice that 

f(vy kl tr( Zl h aMA) ) x 
x tr(r l / is;(A) {/ i ; i ; A) , ^/V( V ,e)}) « -2d6 / f(x){A l (x), yjv{x,e)} A 

where / is any continuous function, v = v(A) and ~ refers to the expansion parameter 
5 of the left hand side in the parametrization of s(e) = v + Ss(e) a (0) + o(<5 2 ), meaning 
that the error is of order 5. We now synchronize^] e ~ 5 and it follows that we can write 
a regulated operator 



x tr([ri/i S;(A )[/i Si(A) , VK])tr(2j^a Jfc (A)) x 

x Xe (t;(A),t;(A0)6^tr([r i ^ (A , ) [/ i ; 1 (A , ) , v^]) x 

x trdAwCA')) (3-10) 

where we could again drop the e dependence in the argument of the volume operator. The 
negative sign in ( |3.10|) stems from the (— i) 2 coming from replacing the Poisson brackets 
by commutators times l/(ih). 

So far everything was true for an arbitrary triangulation. We now apply the operator 
( 3.10| ) to a function / cylindrical with respect to a graph 7 and adapt the triangulation 



to the graph in exactly the same way as in [11, [L2| and as indicated above. Let E(v 



n(v)(n(v) — l)(n(v) — 2)/6 where n(v) is the valence of the vertex v. As we evaluate the 
operator we find out that only those terahedra A in ( |3.1C|) contribute whose basepoint 



v (A) coincides with a vertex v of the graph due to the presence of the volume operators 
in ( 3-10| )- This mechanism is explained in more detail in fTI , I2[ . Moreover, as we take 



e sufficiently small we see that only pairs of tetrahedra contribute which have the same 
basepoint f(A) = v(A'). Combining both observations, we find that we need to sum 

2 This is justified because the classical limit schemes lim e ^o lim<5_>o arL d (lim$_ >o)|e=e((5),e(o)=o are iden- 
tical, both give back the classical magnetic Hamiltonian constraint. We will here neither give the details 
of the dependence e(S) nor a proof but refer the reader to p^j for more details. We will take advantage 
of this fact without mentioning also in later subsections of this paper. 
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only over vertices of the graph and for each vertex v over those tetrahedra A such that 
v (A) = v. Then we find 



Z«Q^lzaj t p ugV(7) v(A)=v(A')=v 



mnp 



jM tr([r i /i Si(A) [/i Si J A) , VKDtrd^.^A)) x 

te{[Tih Spi A')[K^ AI) , V // ^])tr(r / /l Qmri (A'))/7- ( 3 - n ) 



The label "T" on the operator in the first line of ( |3.11| ) is to indicate its dependence on the 
triangulation fjTT], |T2[ which expresses itself partly in the huge freedom of how to choose 
the loops (Xij. This arbitrariness is somewhat reduced in the diffeomorphism invariant 
context that we are interested in because then it does not matter how "large" the loops 
(Xij are as long as the prescription how to attach them is diffeomorphism covariant. See 



TT| , p!2| for further discussion of this point. 



This furnishes the regularization of the Yang-Mills Hamiltonian constraint. 
3.2 Fermionic sector 

In this section we will only focus on the first term displayed in the expression for H]j irac 
in ( |2.2| ). The other two terms can be quantized similarily, for the quantization of K\ we 
adopt a procedure identical to the one used for the quantization of Einstein contribution 



to the Hamiltonian constraint in [n], |T^]. This point is also outlined in appendix A. 

We begin by rewriting the classical constraint using that Ef = \e abc eijke 3 b e^.. We find 
by an already familiar procedure that 

HDirac(N) 

= < f^ M ^ e ^WM- ™ ^» - c.c. 

ZK j Jdet(q)(x) 

(3.12; 

where 6 is an arbitrarily small bu finite parameter. The minus sign comes from moving 
the classical momentum variable to the right as compared to ( |2.2| ). 

The first task is to rewrite ( |3.12 ) in terms of the quantities 9. To that end let /" be a real 



valued, &&su(2) transforming vector field and consider the discrete sum (we abbreviate A/j, 
etc. as / etc.) 

Y,fi(*)frVSi(x)h{x) . (3.13) 



Recall from |7J the definition 9j(x) := / d 3 y^J5(x,y)^i(y) : = lim e ^ 9} (x) where 9\(x 

'Ax, 1 



fd'y^Uy) and 

Xe{x,y) denotes the characteristic function of a box with Lebesgue 



measure e 3 and centre x. We define (d a 9i)(x) := lim e ^ d x a9 e (x) and find 
d*PAx) = Jd 3 y d ^ y) Uy) 

- j d 6 y j- ^(y) = J d J y d y a^(y) 
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since Xe( x ,y) — Xe{y,x) and there was no boundary term dropped in the integration by 
parts because Xe is of compact support. Let us partition £ by a countable number of 
boxes B n of Lebesgue measure e 3 and center x n as in JJJ and interprete ( |3.13| ) as the 
e — > limit of 

J2f?(xn)(T i T> a 9 e ) I (x n )e'i(x n ) . (3.14) 

n 

Substituting for 9 e in terms of £, Q3.14Q becomes 

/ d 3 x J d 3 y f£ f?(x n ) Xn) ^ Xn) } l(ndM(x) + MxnMx)hMy) • (3.15) 

We have not written the Christoffel connection in p.!5| which is needed due to the density 
weight of £ because it drops out in the anti-symmetric sum i[(.) — (.)*] = — c.c] of 
( P-12| ). Now, as e — > (the partition of £ becomes finer and finer) we can replace Xe(x, x n ) 
by 8(x,x n ) and Xe(y,x n ) by 8 Xn , y and fl3.15|) becomes, upon performing the x— integral 
and the sum over x n , 

" d 3 xf-(x)(r t V a ^)(x)Uy) (3-16) 



which is precisely ( ft.12 ) with the proper interpretation of jf . Expression ( |3.16|) is written 
in a form that is well defined on the kinematical Hilbert space which consists of functions 
of 9 rather than £. 

Now, in quantizing expression (|3.13| ) we keep the fermionic momenta to the right and 
replace 9a^(x) by hd/d9A^ which is the proper quantization rule for the 9 variables jTJ. 

Also, we multiply nominator and dominator by 5 3 and replace 6 3 ^Jdet(q)(x) by V(x, S) 
in the denominator which by the standard trick we can absorb into the Poisson bracket. 
Finally we replace the Poisson bracket by a commutator times l/(ih). Labelling the 
regulated operator with the parameter 5, we find on a function / 7 cylindrical with respect 
to a graph 7 with fermionic insertions 9ah at the vertices v £ V(j) 

d 



x 5 3 L4*(x), yJv(x,6)][Ai(x), sjv{x, 5)][(T k V c 9) Alx (v)-—^5 x , v + h.c.]f T (3.17) 

Notice that the sum over all x £ £ already collapses to a sum over the vertices of 7. Next 
we triangulate S in adaption to 7. We have the expansion H s (0, S)9(s(6)) — 9(s(0)) = 
5s a (0)(V a 9)(s(0)). Therefore we just introduce as in the sections before a holonomy at 
various places to absorb the factor of 5 3 and replace V(v, 5) by V v . Thus, 

HiraeiN) = £ N v X 

v(A)=v 



x tr^-^^A)^ 1 ^, VK]) x 



x [(r k [H Sp(A) 9(s p (A)(6)) - 9(v)} A ,^-^ + h.c] 
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■S E *v E 

P ueV(7) u(A)=u 



X tr( T ^ m (A)[^(A)> VK]) x 

x tr^-^^^A), V^])[(n(s p (A)) - y fc («) + /i.e.] 

= : -^Dimc (3.18) 

where the label T reminds us of the triangulation dependence (we have naturally chosen 
the value of 5 in such a way that a) e(S) coincides with the endpoint of the segment of 
e starting at v = e(0) and b) is part of the definition of the triangulation adapted to 7). 
We have defined 

Yi(e) := tT( n H e ae(l))^-rr) and Y&v) := Y(e = v) 

and e : [0, 1] — > £ is a suitable parametrization of the edge e. 

The hermitian conjugation operation "/i.e." involved in ( 3.18|) is meant with respect 



to the inner product on the Hilbert space and with respect to the operator of which the 
first term in Q3.18| ) is the projection on the cylindrical subspace labelled by th graph 7. 
We will return to this issue in the next section. 

Notice that the classical fermionic Hamiltonian constraint is a density of weight one 
and that the operator defined by (|3.18|) precisely respects this because the 9 are scalar 
valued and not density-valued. If we were dealing with the £ instead of the 8 we were 
running into conflict with diffeomorphism covariance at this point. 



3.3 Higgs Sector 

We finally come to regularize the Higgs sector. Especially for this sector a general scheme 
will become evident of how to systematically take advantage of the factor ordering ambi- 
guity in order to arrive at a densely defined operator. 

The term in ( |2.2| ) proportional to (p 1 ) 2 looks hopelessly divergent : even if we would 
manage to replace the denominator by the volume operator we end up with a singular, 
not densely defined operator because the volume operator has a huge kernel. We need a 
new trick as follows : 



We insert the number 1 = [det(e^)] 2 /[ydet(g)] 2 (one) into the kinetic term which ap- 
parently makes the singularity even worse. However, consider the following regulated 
four-fold point-splitting of the kinematical term 

^ H inns, kin (-^0 



1 r a*xN{x) P \x) f Syp\y) f d\{ f*}<> ){u) j £v{ ^ML- )(v) x 



2k 



X Xe(x,y)Xe(u,x)Xe(v,y) 

1 (-2) 2 



o9xN(x)p I (x) i o9yp\y) x 



2k(3\) 2 k 6 

tr({A(u), yVKe)} A {A(u), vV( M ,e)} A {A(u), y/V(u,e)}) X 
tv({A(v), ^V{v,e)} A {A(v), s/V(v,e)} A {A(v), y/V(v,e)}) x 
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X Xe(x,V)Xe(u,x)Xe(v,y) 



(3.19) 



Recall that / d 3 xdet(e l a ) = 4 J e^e 2 A e? A e k = — | / tr(e A e A e) in order to see this. 
Notice that the sign factor in the identity (|3.1| ) has dropped out. We could also have 



used 1 = sgn(det(e*)) det(e Q i)/ydet(^) but then the resulting expression would be less 
symmetric, it is a choice of factor ordering. 

Now we replace p 1 by —ih(K)8 / 8(f) 1 , replace the volume by its operator version and Pois- 
son brackets by commutators times l/(ih) and find, when applying the operator to a 
cylindrical function / 7 , that 

H Higgs,kinV^) fl 

= ^ islv E N(v)Xi(v)Xi(v>) Xe (v,v>)x 
u,«'ev(7) 

x I tv([A(x), Jv(x,e)] A [A(x), y/v(x,e)] A [A(x), y/v{x,e))) x 

x | tr(L4(y), V^e)] A y/v(y,e)] A V^)])te(^«)x £ (!/^') • 

(3.20) 

where X T (v) := |[X^,(f/(f)) + X£(£/(t>))] is the symmetric sum of right and left invariant 
vector fields at U(v) G G. The appearance of X\v) relies on the following consideration, 
explained in more detail in : Instead of p I (x) we consider the integrated quantity p I (B) 
where B is a compact region in S. Now the functional derivative of U(v) with respect 
to (f>i(x) turns out to be meaningless without a regularization of U(v) as well. In |7j 
we use a regularization which takes the interpretation of U(v) as the nontrivial limit of a 
holonomy as e shrinks to e serious. Now the functional derivative can be unambiguously 
performed and since the functional derivative of a holonomy of a connection along an edge 
gives rise to right or left invariant vector fields respectively at the beginning or end of 
the edge respectively it is not surprising that as we remove the regulator on U(v) that we 
obtain a symmetric sum of right and left invariant vector fields. The result is that then 
f{B)U{v) = -ihKXB^X^vpiv). 

Certainly we are now going to triangulate S in adaption to 7 in an already familiar 
fashion and write 



tr([A(x), yjv(x,e)] A [A(x), yjv(x,e)] A [A(x), yjv(x,e)]) 
U^tr(h Si{A) [h^ A) , ^V(v(A),e)]tr(h Sj(A) [h^ {A) , \Jv(v(A), e)] x 



A 



x tr(h Sk{A) [h; k \ A) ,y/V(v(A),e)}) (3.21) 
which results in 

Tft I 

- " E N(p)X I (p)X I (q) X e(p,q)x 



18B 36 

V p,q,r,s£V (7) 



v ' v(A)=r y ' »(A)=s 



X 
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x tr(h Si(A) [h^ {A) , \Jv(v(A), e)]tr(/i Sj (A)[/i~ 1 (A)' \/v(v(A),e)] x 
x tr(/i Sfc ( A )[/i Sfc 1 (A) , \JW(v(A), e)]) 

x tr(/i aj( A0[^(A0> V / ^(^( A, ) ! e)]tr(/i Sm(A0 [/i; m 1 (A , ) , \Jv{v{A'), e)] x 
x tr(/i Sn(A0 [^ A , ))V /5>(A^])/ 7 (3.22) 
since only tetrahedra based at vertices of 7 contribute in the sum J* s = J] a Ja- 

Now we just take e to zero, realize that only terms with v = p = q = r = s contribute 
and find that 

H Higgs ,Um-^ E N^X^X'iv)-^ £ x 

y S i>ey(7) ^ ) v(A)=v(A')=v 



x e i ^tr(/i Si(A) [/i Si J A) , VK,]tr(/i 8j( A)[/i a / (A) , \JV v }tv{h SkiA )[h s ^ A) , \JV V \) x 
x e' mn tr(/i Si(A0 [/ i -J A , ) , V^]tr(/i Sm( AoK4A'). V^M^-W^a')' V^D/" 



7 • 
(3.23) 



The operator ( |3.23|) is certainly quite complicated but it is densely defined ! 



Next we turn to the term containing the derivatives of the scalar field. We write 

E a K b pip k 



q ab Jdet(q) = , \ 1 and E? = e acd e ijk - c d 



and regulate (again we could have chosen to replace only one of the Ef by the term 
quadratic in e l a and still would arrive at a well-defined result at the price of losing sym- 
metry of the expression) 

= i- [ d* X I d*yN( X ) Xe ( X , y ) e «V"V*° ( g ) ^ 



^_(^) 4 / N(x)e ijk V<f> I (x) A {^'(x),\/(x,e) 3/4 } A {A fe (x), e) 3/4 } x 



2k 

-(- 

x I Xe (ar, 2 /)e inm X>^(i/) A {A m (x), V( y , e) 3 / 4 } A {A"(y), V(y, e) 3 / 4 } . (3.24) 

It is clear where we are driving at. We replace Poisson brackets by commutators times 1/ih 
and V by its operator version. Furthermore we introduce the already familiar triangulation 
of £ and have, using that with v = s(0) for some path s 

Ad(k(0,<ft))[E/(s(<ft))] - U(v) = /^0,<ft)[/(s(<ft))/i s (0,<ft)- 1 - U(v) 
= exp(h s (0, 5t)<f>(s(5t))h,(0, dty 1 ) - U{y) 

= exp([l + <fts a (0)AJ [</>(?;) + 5ts a {0)d a <f){v)][l - <fts a (0)AJ + o((5t) 2 )) - U{v) 
= exp(5ts a (0)(d a <p(v) + [Aa, (t>(v))) + o((5t) 2 )) - U(v) = 5ts a (0)V a <f)(v) + o((5t) 2 ), 

(3.25) 
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and with tr(rjTj) = —8ij/2,tT(r I Tj) = —d5u, d the dimension of the fundamental repre- 
sentation of G that 

6 / £>0/(a;) A {A j (x), V(x, e) 3//4 } A {A k (x), V(x, e) 3//4 } 



4 



iniip 



tr(r 7 [Ad(^ m(A) )[f/( Sm (A))] - C/(^(A))]) x 



x to( T jh Sn (A){Krt(A)i V( V (A), e) 3/A })ti(T k h Sp( ^ A - ) {h^ A - ) ,V(v(A), e) 3/4 }) . (3.26) 
Then we find on a cylindrical function 

Hh l99 s4er(N)f, = ^{\)\^) 2 £ A^)x^yK^ m x 

^ t;,v'6V(7) 

x E ^r^tr(r / [Ad(4 Sn(A) )[f/( S „(A))] - 17(«(A))]) X 

v(A)=v v / 

x tr(rA P (A)[< (A) ,^^ x 

x E ■^^'*tr(&[Ad(^ (A0 )[D r ( ar (A'))] - t/MA'))D x 

u(A')=u' ^ > 

x tr^^^)^,), V; 3/4 ])tr(r m / ist(A , ) [/ i ; t ; A/) , V; 3 / 4 ])/, (3.27) 

since only tetrahedra with vertices as basepoints contribute. Thus we find in the limit 
e -> 

HmsssMm^^^ E A^K J V m x 
p ueV(7) 

E ^2^^(r,[Ad(^ (A) )[[/( Sn (A))] - U(v)}) x 

i)(A)=t)(A')=f V ' 

x tr^CA)^^^ x 
x tr(r J [Ad(/i Sr(A , ) )[C/( Sr (A / ))] - U(v)\) x 

x tr(T^ A0 [/i;/ (A , } ^ . (3.28) 

Again, despite its complicated appearence, Q3.28Q defines a densely defined operator. Fi- 
nally the potential term, like the cosmological constant term are trivial to quantize. No- 
tice that certain functions of 4>i(v)(pi(v) can be recovered from polynomials of the func- 
tions [tr(U(v ) n )] m where m,n are non-negative integers. For instance for SU(2) we have 



2 cos(y cj)i(v ) 2 ) = tv(U(v)). Thus we may define for instance a mass term through 

AM' := [arcos(^r 

where the arcos-function is for the principal branch and is well-defined because the ar- 
gument takes values in [—1, 1] only. Thus, by this rule all polynomials in 4>i(v) 2 become 
actually bounded functions of U(v). This is not an unknown phenomenon, the same hap- 
pens when one replaces the Yang-Mills action by its regularized Wilson action on a fixed 
lattice (our lattice, the triangulation, is not fixed, it "floats" with the state). Therefore, 
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taking <pi(v) 2 as expressed through those products of traces we find 

H cosmo (N)f^ = -jjL. £ N v V v f y . (3.29) 

This furnishes the quantization of the matter sector. Notice that all Hamiltonians have the 
same structure, namely an operator which carries out a discrete operation on a cylindrical 
function, like adding or subtracting lines, fermions or Higgs fields, multiplied by the Planck 
mass and devided by an appropriate power of the Planck length which compensates the 
power of the Planck length coming from the action of the volume operator. It follows 
that in this sense the matter Hamiltonians are quantized in multipla of the Planck mass 
when we go to the diffeomorphism invariant sector. 



3.4 A general regular izat ion scheme 

In this subsection we describe a recipe by means of which a fairly large calss of Hamiltonian 
densities of weight one which are diffeomorphism covariant and coupled to gravity can 
be turned into densely defined and, as we will see later, anomaly-free operators on the 
Hilbert space that we have defined. The resulting expression does suffer from a factor 
ordering ambiguity but not from a factor ordering singularity. 
The restrictions on the Hamiltonian density are as follows : 

a) The matter canonical momenta P of the theory are scalar densities of weight one 
and matter configuration variables are scalars (they may transform non-trivially under 
SU(2) x G). In case that matter is a priori described by tensors, turn them into internal 
SU(2) tensors by means of the triad and co-triad, the corresponding canonical transfor- 
mation will add to the gravitational connection a piece K l a which is a real valued one-form 
and transforms homogenuously under 577(2); thus the reality of A l a is preserved under 
this canonical transformation. Other cases require a special treatment (for instance the 
case of the fermion fields). 

b) Furthermore, it is assumed that all covariant derivatives are with respect to A l a ,A^, 
act only on configuration scalars Q and are of first order only so that no Christoffel con- 
nection is needed. In case that the covariant derivative is a priori given in terms of the 
spin connection and/or acts on a tensor, write the tensor as before in terms of (co)triads 
and the canonical configuration scalars. If the covariant derivative is D a with respect to 
A a we just use that V a e\ = e^ k K 3 h e k . If it is D a with respect to the spin connection, use 
that e l a is annihilated by D a and that (V a — D a )vi„j = K k [eikiVi„j + .. + eju v i..i\- We see 
that all covariant derivatives can be cast into the desired form up to underived factors 
of K l a which we write as a Poisson bracket {A l a , {H E (1), V}} where H E is the Euclidean 



Hamiltonian constraint [0, 112 



Restrictions a),b) are just in order to state the theorem below in a compact form. The 
case of higher derivatives (as they actually occur in the Riemann curvature) just require 
a case by case analysis. 

Consider then a general Hamiltonian density which is local and consists of monomials 
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of the form (we suppress all indices and contractions other than tensor contractions) 



H m>n (x) = [P(x)] n E a ^x)..E a -(x)f m , n [Q] ai .. am (x) r ) -. (3.30) 

/det(q)(x)} n+m - 1 

Here f m>n [Q] is a tensor of density weight zero which is independent of e l a ,P and is a 
polynomial consisting of sums of terms involving covariant derivatives of Q of first order 
and underived factors of K l a or F l ah such that the total number of their covariant indices 
is m. The denominator accounts for the fact that the Hamiltonian is a density of weight 
one. 

Expression ( |3.30| ) defines the most general basic building block of the Hamiltonians under 
consideration, that is, every Hamiltonian in the class that we have defined is a linear 
combination of these. 



Theorem 3.1 (Structure Theorem) Any Hamiltonian constraint of the form ( 3. 3Q) 



can be turned into a densely defined operator on 7i which is diffeomorphism covariantly 
defined and anomaly free. 

Proof : 

In case m + n = we are done because upon quantization we just need to triangulate S 



and replace ydet(q)(x) by V x . 

Consider then the case m + n > 0. In order to regulate (|3.30|) we will need m + n — 1 point 
splittings for the n momenta P, E. This will require m + n — 1 regulated ^-distributions 
X e /e 3 . Each factor factor 1/e 3 can be absorbed by replacing l/waet(g)(x) by l/V(x,e) 



but then we cannot simply replace this by its operator version because V(x, e) is in the 
denominator. 



Now multiply ( |3.30|) by the number 1 = [| det(e^)|/wdet(g)] fc and introduce k > more 



point splittings. We have a power of 3k co-triads e l a in the nominator and a power of 



n + m + k — 1 factors of e 3 ydet(g) in the denominator. We replace each e 3 y det(q)(x) by 
V(x, e) following the standard trick. 

Now let e := sgn(det(e^)) = (e) 3 . We have the classical identity, using (|3.1| ) 



det(<)(*)| = ^e{x)e abc e l]k {e\e{e k c ){x) = ^e ab % jk [(eel)(eei)(ee k c )}(x) 

' ^^{A^x), V(x, e)}{Al(x), V(x, e)}{A k c (x), V(x, c)} (3.31) 



and therefore each e l a is worth a factor of V(x, e) in the nominator within a Poisson 
bracket. Now choose k large enough until 3k>n + m + k — 1, i.e. 2k > n + m — 1. By 
suitably point splitting the various factors we get 3k factors of the form 

{A\{x\ V(x, c)} {Al(x), V(x, e) 1 -^^ 1 } 



V(x,e) 



n+m+k-i n+m+k—1 



3k 



3k 



and thus have managed to produce a net positive power of volume functionals for the 
point-split classical Hamiltonian density in each of the 3k Poisson brackets. We now 
choose the arguments x, y of the various Xe(%,y) so that the limit e — >• gives a non- 
vanishing result only if all arguments coincide. We triangulate £ in the fashion outlined 
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in the previous subsections and replace A in J A tr({A, V} A {A, V} A {A, V}) by the 
holonomies along the edges of the triangulation. Finally we replace P, E, ordered to the 
right, by the corresponding functional derivatives, Poisson brackets by commutators and 
the volume functional by its operator version. The result when applied at finite e to a 
function cylindrical with respect to a graph only gives contributions at an m + n + 3k 
tupel of vertices (or edges) of the graph and when sending e to zero all vertices of the 
tupel have to coincide in order to give a no n- vanishing result. This shows that we find a 
densely defined operator. 

To see that it is anomaly- free we just note that the resulting operator is of the type 
to which the theorem of the section on anomaly-freeness applies. 
□ 

Remarks : 

1) We note that the density weight of one was crucial (besides the fact that the 
integrated operator is only diffeomorphism invariant if the density weight is one) : 

If it would have been higher than one then we needed n + m + k — 1 point splittings but 
in the denominator we have a power of y / det(g) which is smaller than n + m + k — 1 and 
therefore even the regulated operator blows up at least as 1/e 3 . If it was less than one 
then by a similar argument the regulated operator vanishes at least as e 3 which is trivially 
always zero. 

2) The proof shows precisely the sources of the factor ordering ambiguity : 

a) That we chose the momenta to the right was essentially forced on us because we want 
to obtain a densely defined operator : if the functional derivatives act on f m ,n(Q) then 
in general it will not be true any longer that the operator only acts at the vertices of 
the state but at all vertices of the triangulation which are infinite in number and so the 
resulting state would not be normalizable. 

b) As long as 2k > n + m — 1 we can have arbitrarily large k and still get a well-defined 
result. Surely, minimal k is desirable to obtain a simple result. 

c) We could have absorbed different powers of V(x, e) into the various Poisson brackets, 
however, all powers must add up to n + m + A; — 1. 

d) We are free to take advantage of the classical identity 

E?(x) = \e abc e l]k e{{x)e k c {x) = ^e ahc e ijk {^{x), V(x, e)}{A c {x)\ V(x, e)} (3.32) 

to lower the necessary value of k if desirable since each of the m factors of E is worth a 
power of two of V(x, e). Of course, it may be true that the value of k must be at least one 
in order to get a densely defined operator (if k > then for sure the resulting operator 
will act only at vertices of the graph as we proved). 

3) The theorem works the same way in any dimension d > 2 because the critical con- 
dition dk>n + m + k — 1 can be satisfied by some k for any value of n, m only for d > 1. 
However, in one spatial dimension all tensors are densities and in zero spatial dimension 
we do not have a field theory so that the theorem does not apply in these cases anyway. 
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4 Consistency 



In this section we will perform the required consistency checks necessary to show that we 
really constructed covariantly defined, anomaly- free, linear operators through their action 
on cylindrical functions which is non-trivial in the sense that it has a non- vanishing kernel. 



4.1 Cylindrical Consistency 

Notice that, just like the gravitational Hamiltonian constraint, the matter Hamiltonian 
constraints are actually defined as a family of operators {Hi} where I is a compound 
label consisting of the graph, colours and spins of its edges and fermionic and Higgs 
representations of its vertices, that is, I = (7, [j, n,p\, [c, C, q\) =: (7, A). The set of labels 
/ is an uncountably infinite one because the set of piecewise analytical graphs of E has 
this cardinality. Still this set allows for a nice and controllable orthogonal decomposition 
of the Hilbert space. In analogy with the source-free CclSC, clS the reader can easily prove 
himself given the measures defined in section 2, we have : 

where the first direct sum is an uncountable one while the second is a countable one, 
running over the possible colourings of the graph with the various compatible irreducible 
representations, compatible in the sense that there exist projectors which render the 
associated cylindrical functions into gauge-invariant ones. The Hilbert space 7i 1 is infinite 
dimensional and is the completion of the space of functions built from spin-colour-network 
states on 7 while % is a finite dimensional vector space with dimension equal to the 
number of linearly independent projectors on gauge invariant functions compatible with 
the colouring A of 7. 

As the decomposition in ( flip is direct there exist orthogonal projections P^xH.*-^ Ti^ 



Cylindrical consistency now means that the family of operators {H ^} is a family of 

projections of a single operator H defined on 7i such that Hp^ % = H^ The necessary 

and sufficient condition for this to be the case is that H , % = whenever 7 7^ 7' or 

A 7^ A'. But this is the case by construction if we simply define H := XL \H^ This 
suffices to prove consistency. 



4.2 Diffeomorphism-Covariance, Continuum Limit 
Self-Adjointness and Positive Semi-Definiteness 

We begin with diffeomorphism covariance of the family of operators obtained and the 
final continuum limit given by the infinite refinement of the triangulation. 
Let $ be the topological vector space constructed in || [Tj of finite linear combinations of 
spin-colour-network functions. By $' we mean the continuous linear functionals on $ and 
denote by Hg(N) the regulated Hamiltonian constraint where the parameter 5 > is a 
regularization parameter expressing the fact that the loops attached are finite in size, that 
is, we did not take the continuum limit yet. Let / = J2i C1T1 G $ be a cylindrical function 
where Tj are the spin-colour- network states. As shown in || 0, the following object makes 
sense as a distribution on $ : [Tj] := J^TeiT^T, where {Tj} := {U((p)Tj : G Diff(S)} 
is the orbit of the vector Tj under the diffeomorphism group. Moreover, as one can 
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show, the [Tf] provide an orthonormal basis on the diffeomorphism invariant Hilbert 
space, that is, any \l/ as above is a linear combination of those. Furthermore, we have 
[/] := J2i c i[Ti], that is, every / G $ gets averaged term- wise after decomposing it into 
spin-colour-network states (the reason for this term-wise averaging is explained in || |7j , 
we also neglected here some technical details which one can also find in those papers). 
The definition of the inner product on the space of diffeomorphism invariant distributions 
is given by < /, g >Diff'-— [/](<?) where the latter expression means the evaluation of the 
distribution [/] on the test function g [§, [/]]. 

A diffeomorphism invariant distribution \I/ G $' is a solution of the Hamiltonian constraint 
provided that 

*{H s {N)f) = (4.2) 

for any / G $ and iV G S (the space of test functions of rapid decrease). The striking 
feature of (|4.2|) is that it is independent of the value of 5 ! The underlying reason for this 
is the diffeomorphism covariance of the graph- dependent triangulation prescription. This 
is proved in Jl2|, ^] for the gravitational part of the Hamiltonian constraint and the same 
reasoning applies to the matter coupled case as well and will not be repeated here. 
The limit 5 — > is therefore already performed in (|4.2|) (see also [0)- Notice that, just 
as in the continuum limit of a regulated Euclidean path integral of constructive quantum 
field theory, we take the limit after integrating. Indeed, the limit before integrating does 
not exist in the L 2 sense (for the same reason that the generator of the diffeomorphism 
constraint does not exist 0). 

Next we turn to the following issue : 

If we set iV = 1 and the spatial metric has signature (+, +, +) then the matter contribu- 
tion to the integrated Hamiltonian constraint, at least for the Yang-Mills and the Higgs 
part, is classically a manifestly non-negative and diffeomorphism invariant functional. 
Since upon replacing the dynamical gravitational field by some classical background field 
this functional plays the role of the matter Hamiltonian in that background field, a natural 
question to ask is whether it should not be promoted to a positive-semi-definite, diffeomor- 
phism invariant operator in the quantum theory (with dynamical quantum gravitational 
field). In fact, this question is even more natural to ask in view of some kind of "quantum 
dominant energy condition" on the Hamiltonian matter density (which equals the Hamil- 
tonian constraint) because it is the energy density component of the energy momentum 
tensor (see |IU| for a first attempt towards a quantum formulation of a dominant energy 
condition in the quantum theory). 

While imposing positivity is then very natural from this point of view it is, in fact, unnat- 
ural from another point of view : namely, if the matter Hamiltonian operator is positive 
semi-definite, then it is hard to imagine how that can be true if not the Hamiltonian den- 
sity, when integrated over any compact region of S, also becomes a positive semi-definite 
operator. In particular, the matter Hamiltonian constraint when integrated against a 
non-negative lapse function should also be a postive semi-definite operator. However, 
then this latter operator will be automatically symmetric (it even would have self-adjoint 
extensions, at least its Friedichs extension). Now, by arguments explained in [pql , a sym- 
metric Hamiltonian constraint operator is in danger of being in conflict with the task of 
faithfully implementing the Dirac constraint algebra. Most certainly, one expects a quan- 
tum anomaly in this case. Therefore, if one wants an anomaly-free quantum constraint 
algebra then it seems that one should not insist on a positive semi-definite Hamiltonian 
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constraint operator. One might think that this can be accomplished, while keeping the 
matter Hamiltonian constraint positive, by having a non-symmetric gravitational Hamil- 
tonian constraint but, at least on solutions, the gravitational and mater contribution just 
equal each other up to a sign and therefore necessarily the gravitational piece is also 
symmetric if the matter piece is. 

We will leave the resolution of this puzzle for future investigations and just mention 
that : 

1) As we will see in the next section, the Dirac algebra is indeed not faithfully represented 



one could, however, use the arguments of |27| to improve this), but still the algebra is 



non-anomalous in the sense that we obtain a consistent quantum theory, 

2) The full matter Hamiltonian operator, as it stands in this paper, is not symmetric but 
at least it is diffeomorphism invariant, 

3) The electric piece of the Yang-Mills Hamiltonian and the term bi-linear in the Higgs 
momenta are indeed essentially self-adjoint and positive semi-definite operators on the 
diffeomorphism-invariant Hilbert space. 

In the sequel demonstrate 2) and 3). 

Denote by H™ := H™(1) the matter Hamiltonian constraint evaluated at unit lapse. 
We define the diffeomorphism invariant analogue of H™ in the continuum by 

< [T T ], H m [Tj] > Dtff := \xm[T I )(Hg l T J ) (4.3) 

where on the right hand side we have the evaluation of a distribution on a test function. 
Now, since H^(N)Tj G $ we have by definition [Tj)(H^Tj) =< [Tj], [HfTj] > Diff but 
since [Tj] is diffeomorphism invariant, this number, as explained in ( |4.2|) , does not depend 
on 5 > (the size or shape of the loop attached) and so is a constant. The limit is 
therefore trivial and so 

< [Tj], H m [Tj] > Dlff = [Tj](H?Tj) (4.4) 
for any 5 > 0. This displays H m as a diffeomorphism invariant operator. 



We will now demonstrate that ( |4.4|) is not even symmetric : 

We restrict ourselves to the case that we couple only a pure Yang-Mills field. We just 
need one counter-example : Let Tj be a spin-network state based on a graph 7 with 
three edges and two tri-valent vertices only one of which, v , is such that the tangents of 
the edges are linearly independent at it. Colour the edges with suitably high irreducible 
representations of SU(2) x G such that Hy M ^Tj is a sum of spin-colour-network states 
Tji each of which depends on an extended graph 7' of which 7 is a proper subset, that 
is, 7' = 7 U A where A is a tetrahedron at v defined by the triangulation. Now choose 
Tj := Tji for one of the J' so that [Tj'](Hym &Tj) 7^ 0. On the other hand, if we apply 
Hym s to Tji then we get a linear combination of spin-network states which depend on 
graphs each of which includes 7 but is even larger than the graph 7' which underlies J'. It 
follows that [Tj}(Hy m sTj') — thus contradicting symmetry. A similar argument reveals 



that the Dirac and Higgs Hamiltonians as defined by ( |4.3|) cannot be symmetric 



Let us then conclude this subsection with showing that the electric piece of the Yang-Mills 
Hamiltonian is a densely defined, positive definite essentially self-adjoint operator even on 



26 



the auxiliary Hilbert space TC. A similar argument applies to the p T p T part of the Higgs 
Hamiltonian. 

Let Of B := (h s [h~ 1 ,V])AB for any segment s where A, B denote SU{2) indices. Using 
that the adjoints on H of {hs) ab^ V are given by {hs) ab^ V respectively and the unitarity 
of SU{2) we compute 



(Of) T = [V,(h 7 i)c B ](h s ) A c 

= -{{{h s ) BC {h- l )cA, V] - ^BcKK'hA, V]) 

= +Of A - (4-5) 

Next, using that Xj s commutes with Of B and that (Xf)t = —Xj s is anti-self-adjoint on 
H we find for Of*' 1 := dj B X{ that (6 s AB ' / ) t = -Of AJ ■ 

Putting everything together we find that when setting N — 1, H E M =: J2v k E H E Mv with 
k E > that 

= E E^fW^EECTE^] (4.6) 

vGsDs' A,B;I A,B;I vEs v<=s 

whence we have displayed H in the form of a sum of operators of the form A' A, that 
is, it has positive semi-definite and symmetric projections. Since these projections map 
Ti-^x m t° itself (the dependence on the segments s(e) involved in the definition of H E M 
drops out because of gauge invariance, see 0]), it follows that the family of projections 
defines a positive and symmetric operator on 7i which therefore has self-adjoint extensions 
(actually, by methods similar as for the volume, area and length operators |21], ^2], |23], [24)1 



one can show that each projection is essentially self-adjoint and so is the whole operator, 
therfore the extension is unique). Notice also that the electric part does not depend on the 
regulator 5 any longer so that the continuum limit is already taken. It therefore projects 
to a strongly diffeomorphism invariant self-adjoint and positive semi-definite operator on 
the Hilbert space defined by (f4.3| ) by general theorems proved in 0. 
We mention that it is conceivable that the positivity of the electric piece could be sufficient 
to establish positivity of the full Yang-Mills Hamiltonian if the Kato condition applies to 
the symmetrically ordered (but not positive) magnetic piece. 

4.3 Anomaly- Freeness 

We will understand the term "Anomaly-free" in the sequel to mean that we have a con- 
sistent quantum theory, namely, the commutator of two Hamiltonian constraints vanishes 
when evaluated on a diffeomorphism invariant distribution. We do not mean that the 
commutator equals a certain operator that is proportional to a diffeomorphism generator. 
The difficulty in achieving this more ambitious goal is two-fold as pointed out already in 
HI : First, not even the generator of a one parameter family of diffeomorphisms exists 



in the representation that we have chosen since the associated representation of the dif- 
feomorphism group does not act strongly continuously on the Hilbert space. Secondly, 
since {H(M),H(N)} = V a (q ab (M )b N - MN b ) where V a is the classical generator of the 
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diffeomorphism group, one needs to make sense out of an operator that somehow corre- 
sponds to q ab which is not at all obvious to construct. By the structure theorem proved 
in the previous section we know that there is a chance that there exists a well-defined 
operator corresponding to the product q ab Vb. Indeed, in |26| such an operator is con- 



structed and part of it actually does generate diffeomorphisms ! However, the way that 
it results from computing the commutator is rather non-standard and has to do with rea- 
sons deeply rooted in the structure of the Hilbert space TC. Basically, one can show that 
the commutator is weakly equivalent to that operator because both are zero operators on 
diffeomorphism invariant states. 

This is almost equivalent to showing anomaly-freeness in the weak sense as stated 
above with which we content ourselves here and which by itself is also a non-trivial task. 
The computations are very similar to the vacuum case so that we refrain from displaying 
all the details. The interested reader is referred to |l2j to fill the gaps. 
To begin with, notice that when evaluated on a cylindrical function, the Hamiltonian 
constraint of both gravity and matter is a sum of terms of the structure N V H V where H v 
is an operator built from 

1) holonomies of segments of the underlying graph which start at the vertex v 

2) gravitational Volume operators which act only on holonomies along segments starting 
at v 

3) Yang-Mills Laplacian operators which act only on holonomies along segments starting 
at v 

4) Fermion field and Higgs field derivatives which act only on fields located at v 

5) Fermion field and Higgs field insertions at the vertices of 7. 

The crucial point is that all the terms involved in H v involve a factor of the form V v or 
more generally V™, n > where again this notation means the volume operator for an 
arbitrarily small neighbourhood of the vertex v. Notice that if / is a function cylindrical 
with respect to a graph 7 then H(N)f = H 7 (N)f = J2veVM N v H v>1 f is in general 
a function cylindrical with respect to a graph 7(f) which contains additional vertices, 
but these vertices are co-planar (arising from loop-insertions due to the gravitational 
or magnetic part of the Yang-Mills Hamiltonian constraint) or co-linear (arising from 
fermion field insertions due to the Dirac Hamiltonain coinstraint). Therefore, V v > = for 
v' G 7(f) — jiv G V(j). We therefore conclude that just like in the source free case 

[H(M),H(N)]f 
= ( M vN v - M v N v ,)H l{v) yH^ v f 

v,v'eV(y) 

= 2 ^ {M V 'N V — M v N v >)(H 7 ( v ) )V rHy jV — H 7 ( v rj jV H %v t)f 

v,v'(zV(y) 

1 . . . . 

= - ^ (M V >N V — M v N V ')(H 7 ^yH JjV — H 1)V /H^^ v )f 

Vy^V' ,v ,v' &V ("f) 

(4.7) 

In the last step we used that the terms with v = v ' trivially vanish while for the terms 
with d/d' the local character of the operator H ltV , performing changes of 7 only in a 
neighbourhood of v makes it commute with H^r^y. It is then easy to see, by the same 
argument as in [f[2]| , that the two functions Hj^yH J}V f, H^yH^^^f are related by an 
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analyticity preserving diffeomorphism. 

This furnishes the proof of anomaly-freeness. 



4.4 Solutions of the Diffeomorphism and Hamiltonian constraint 

The general solution of the Diffeomorphism constraint for theories including Fermions and 
Higgs fields is constructed in 0. As for the pure gauge field case ||, they are elements 
of$'. 

We can now look for solutions to the Hamiltonian constraint. Recall (|12||) that a 
solution to the Hamiltonian and diffeomorphism constraint is a distribution \P G that 
is, a continuous linear functional on the space $, the finite linear combinations of spin- 
colour-network states such that 

1) *[U(<p)f\ = *[/] V ip G Diff(E), / G $ and 

2) V[H(N)f] = V N G 5(E), / G 

Here, as usual, 5 is the Schwartz space of functions of rapid decrease. 
The complete set of solutions can now be precisely characterized as follows : 
Let R := UAr e ,sRan(i^(iV)) be the union of the ranges of H{N) on <£> and let S := R 1 - be its 
orthogonal complement in $. Next, for each s G S there is a decomposition s = J2i s iTi 
of s into spin-colour-network states. Consider for each Tj its orbit {T/} := {U(ip)Tf, <p G 
Diff(E)} under diffeomorphisms and construct the distribution [Tj] := J2t&{T;}T G 
That this is still an element of $' follows from the fact that above we took the orthogonal 
complement in <3> and not in 7i. Define [s] := J2i s i[Ti]- Then the the complete space 
of solutions to both constraints is given by V p h ys the (infinite) linear combinations of ele- 
ments of the set {[s] s G S} C 

We see that we know the space of solutions once we know S. A compact algorithm to 
describe S as for the source-free case is not available at the moment for the matter coupled 
case so we restrict ourselves to displaying trival solutions which is enough to show that 
the theory is not empty : 

Trivial solutions are, for instance, distributions which are defined by graphs all of whose 
vertices are co-planar, just because the volume operator annihilates such states so that 
they cannot be in the image of the Hamiltonian constraint. The dependence on Fermion 
and Higgs Fields of these solutions is completely arbitrary. 

While we do not have an explicit algorithm for the construction of the kernel at the mo- 
ment, it is clear that it exists and therefore one can construct a physical inner product 



and strict Dirac observables as outlined in [26]. 
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A Canonical real-connection formulation of Einstein- 
Dirac Theory 

In this section we wish to derive the canonical action principle for general relativity coupled 
to Dirac fields in manifestly real form and in terms of the canonical pair (A a = T a +K a , E?) 
introduced in section 2. This has not been done so far in the literature because either 
one was interested in the standard Palatini formulation or in the Ashtekar formulation, 
the latter involving complex-valued connections which are difficult to deal with in the 
quantization programme. The progress that one is able to make with the real-valued 



variables for the source-free case |12| motivate to derive a similar form of the action in 
the matter-coupled case. 

Let us begin with the (massfree) Dirac action in covariant form 

S Dir ac =\ \ dV-det(<?)[* 7 XV a * - W^^Jf] (A.l) 

2 JM v 

where 7° are the Minkowski space Dirac matrices, e a a are the tetrad fields and \1/ = (ip,T}) 
is a Dirac bi-spinor and \I/ = (\I/*) T 7 its conjugate. Here ip = (ip A ) and r\ = (fjA 1 ) 
transform according to the fundamental representations of SL(2, (D) and are scalars of 
density weight zero. The covariant derivative V a is defined to annihilate the tetrad e^, 
that is, we are using the second order formalism. 

In order to put ( |A.1| ) into canonical form we take M = R x E, let T a be the time foliation 
vector field of M and denote by n a the normal vector field of the time slices S. Then the 
tetrad can be written e" = e° — n a n a with e"n„ = e a n n a = so that e" is a triad and 
rj al3 n a i>f3 = —1 is an internal unit timelike vector which we may choose to be n a = — 5 a ,o 
(rj = diag(— , +, +, +) is the Minkowski metric). Finally, inserting lapse and shift fields 
by {dt) a = T a = Nn a + A^ a with N a n a = one sees that the action can be written, after 
lengthy computations, in terms of Weyl spinors as (using the Weyl representation for the 
Dirac matrices, for instance, to expand out various terms) 



% r r 1 T a — N a 

Soirac = - dt d 3 xN^det(q)[ (^V+V + V* VaV - c.c; 

2 jr Je v JM 

+e^V iV +^ - rfviV'V - c.c.)} (A.2) 

where the c.c. in (* + c.c.) stands for "complex conjugate of *. Here we have defined 
e a a = (0,e") and abused the notation in writing e" = (e* = 0, sgn(det(e))e"), 3 Oi are the 
Pauli matrices, ift := (4>*) T and is the self-dual respectively anti-self-dual part of 
V a in the Weyl representation. More precisely, = d a + uf, = —iujU^, u 3 a + = 

— l/2ejkiu>a l+ , = §0^°^ ~ ^ a/3 -Y5 u a 5 )i ^a' = an d ls the spin-connection 

of e a a . The unfamiliar reader is referrred to the standard literature on the subject (for 
instance p8f). 

It is easy to see that the spatial part of u 3 a + is just given by \A 3 a € where A{ € = T{ + iK{ 

is the complex-valued Ashtekar connection. Denoting T>®ip = (d a + Aj~Tj)ij) and T> a rj = 

(d a + Aj € Tj)r] with Tj = (Pauli matrices) and A{ € = T a u a + , ip = T a d a ip we end 

up with 

S Dirac = I J dt J d 3 x ^det (q) [(ip^ip + rfr] - c.c.) 
3 the sign factor is a possible choice because it does not appear in q a (, = e l a e\ 
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-(-(Af^Tjip + AfrfTj-rj - c.c.) + N a (^Vfifj + rfV^v - c.c.) 
+iVsgn((det(e)K(-VV i P£V + V^i^aV ~ c.c.))} . (A.3) 

Let us now introduce D a ip = (d a + TjT :) a )ip, Ef = det(e^)e", A\ := $t{A\ ' ) then we see by 
explicitly evaluating c.c. that 



c.c. 



S Dirac = I dt c/ 3 x^/det(g)[(^V + rfr] 

-{-2A\{^t^ + r^Tj-n) + N a (^D a ij + rfD a ri - c.c.) 

+JV_2=([-^V < D a V + V^iDaV ~ cc] + 2[K a , E a ] j (if)Wjil> - tjV^jjJM) 
'det(g) 



This is the 3+1 split Dirac action that we are going to combine with the 3+1 split Einstein 
action to obtain the desired form in terms of (A l a , Ef). 

We come to the Einstein action. In contrast to [18| we also take the second order 
form of the Palatini action (that is, we let the gravitational connection be the one that 
annihilates the tetrad from the outset). Otherwise we can take over the results from [IS 
and arrive at SEinstein = $1(Se) where is the self-dual part of S Einstein which in our 
notation is written as 

S + E = -= jdtf d 3 x[-zAfE« - (iAfV^E^ - iN%r(F^ b E b ) + N tr(if b [£", E b ])] 
k J J 2Jdet(q) 



(A.5) 

where F c denotes the curvature of A € and k the gravitational coupling constant. Com- 
puting the real part reveals 

SEinstein = \ = j dt j d 3 x[KlE^ - (—A 3 t [K a , E% + 2N a D [a Ki } E b 

A ti((\K a ,K b ]-R ab )[E a ,E b ])] (A.6) 



2 v /det(g) 
where R a b is the curvature of e\. 

Thus, putting both actions together, we find that the gravitational Gauss constraint is 
given by (no other matter contributes to it) 

Gj = ~[K a , E% + iy/detiq^Trf + Tj^rj) . (A.7) 

We can now perform a canonical point transformation on the gravitational phase space 
given by (K l a , Ef) —* (A l a , Ef) (the generator turns out to be / d 3 xT l a Ef as one can 
explicitly check) and we must then express the constraints in terms of A l a . Let us therefore 
introduce the real- valued derivative V a ip := (d a +A 3 a Tj)i/} and denote by F ab the curvature 
of A l a . Using that D a Ef = we can immediately write 

Q j = -V a E] + z^/det^V^ + rfTtf] . (A.8) 

Next, we expand F a b in terms of T a , K a , use the Bianchi identity tr(R a bE b ) = and find 
that the vector constraint V a , the coefficient of N a in Sriirac + SEinstein is given, up to a 
term proportional to Qj, by 



V a = tr(F ab E b ) + -^/det(q)(^V a ^ + rfVj) - c.c). (A.9) 
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Finally, let as in the source-free case 



H E = ^tr(F ab [ ^£l) (A.10) 

which has the interpretation of the source-free Euclidean Hamiltonian constraint. Fur- 
thermore, let 

Ha ■= -H E + £-ix([K a , ^]^=|) (A.11) 
2k Vdet(g) 

which in the source-free case would be the full Lorentzian Hamiltonian constraint. Then 
the Einstein contribution to Hamiltonian constraint of S^irac + SEinstein is given by 

H = H G - ^D a tr([K b , E b ] E " ) =:H G + T. (A.12) 
2k Vdet(g) 

Notice that in the source-free case the correction T of H to H G is proportional to a Gauss 
constraint and therefore would vanish separately on the constraint surface. However, in 
our case, using the Gauss constraint (|A.7|) we find that 

T = --([#«, E a Y - Ep a )Jj (A.13) 

where we have defined the current Jj = ip^tjjip + V VjV- On the other hand, writing also 
the Dirac contribution to the Hamiltonain constraint in terms of T> a rather than D a and 
combining with H we find that the first term on the right hand side of ( |A.13j ) cancels 
against a similar term. We end up with the contribution C from both the Einstein and 
Dirac sector to the Hamiltonian constraint which is given, up to a term proportional to 
the gravitational Gauss constraint, by 

E a A 



C = H G + —=J—{V a (Jdet(q))J j 



+i y Jdet{q)[^c j V a il) - rfoPari - c.c] - AT^det^^V - rfv)) • (A.14) 

In order to arrive at (|A.14j) one has to use the Pauli matrix algebra OjOk = Sjklsu(2)+i^jkiCi 
at several stages when computing c.c. Notice that we can write ( A.14j) also in terms of 
the half-densities £ = ^det (q)^jj, p = ^det(q)r] by absorbing the ^det(g) appropriately 
and using that V a det(q) = 0. We find 

c = h g + E " (^(ev^+pVp) 

2Jdet(g) 



+i[evPai - P^PaP ~ cc] - K> a @£ - pt p )) . (A.15) 



Note also that iJdet(q)[ip^ip ~ ipty] = — so that our change of variables is 
actually a symplectomorphism ! 

This is the form of the constraint that we have been looking for : up to K l a we have 
expressed everything in terms of real- valued quantities and the canonically conjugate 
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pairs (£, (p, ip). Now, let as in the source-free case denote V = j d 3 x^Jdet(q) the total 
volume of £ and H E (1) = J d 3 xH E (x). Then it is still true that K{ = -{A{, {V, H E (1)}} 



and since V,H E (1) admit well-defined quantizations [Tj we conclude that despite its 
complicated appearance ([A. 15 ) admits a well-defined quantization as well. Note that 
if we had not worked with half-densities £, p but with the ip,r] then, while i^Jdet(q)ip 
is the momentum conjugate to ip, the gravitational connection would get a correction 
proportional to ie\[^ijj + rfrj\. Thus we would have had to admit a complex connection 
which would be desasterous as the Hilbert space techniques [3| |], |||, || would not 

be at our disposal. Therefore the strategy of working with ^'s as advertized in || is 
not suitable for quantizing the Einstein-Dirac theory and we are forced to adapt the 
framework developed in @. 
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